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ing structure of a recent theory of the Ising-nematic transition in metals, and discuss 
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1 Introduction 

The past couple of decades have seen vigorous theoretical activity on the quantum 
phases and phase transitions of correlated electron systems in two spatial dimen- 
sions. Much of this work has been motivated by the cuprate superconductors, but 
the list of interesting materials continues to increase unabated QJ. 

Methods from field theory have had a strong impact on much of this work. In- 
deed, they have become part of the standard toolkit of condensed matter physicists. 
In these lectures, I focus on two classes of strong-coupling problems which have 
not yielded accurate solutions via the usual arsenal of field-theoretic methods. I will 
also discuss how the AdS/CFT correspondence, discovered by string theorists, has 
already allowed substantial progress on some of these problems, and offers encour- 
aging prospects for future progress. 

The first class of strong-coupling problems are associated with the real-time, fi- 
nite temperature behavior of strongly interacting quantum systems, especially those 
near quantum critical points. Field-theoretic or numerical methods often allow ac- 
curate determination of the zero temperature or of imaginary time correlations at 
non-zero temperatures. However, these methods usual fail in the real-time domain 
at non-zero temperatures, particularly at times greater than h/ksT, where T is the 
absolute temperature. In systems near quantum critical points the natural scale for 
correlations is Ji/ksT itself, and so lowering the temperature in a numerical study 
does not improve the situation. 

The second class of strong-coupling problems arise near two-dimensional quan- 
tum critical points with fermionic excitations. When the fermions have a massless 
Dirac spectrum, with zero excitation energy at a finite number of points in the 
Brillouin zone, conventional field-theoretic methods do allow significant progress. 
However, in metallic systems, the fermionic excitations have zeros along a line in 
the Brillouin zone (the Fermi surface), allowing a plethora of different low energy 
modes. Metallic quantum critical points play a central role in many experimental 
systems, but the interplay between the critical modes and the Fermi surface has not 
been fully understood (even at zero temperarture). Readers interested only in this 
second class of problems can jump ahead to Section|7] 

These lectures will start with a focus on the first class of strong-coupling prob- 
lems. We will begin in Section|2]by introducing a variety of model systems and their 
quantum critical points; these are motivated by recent experimental and theoretical 
developments. We will use these systems to introduce basic ideas on the finite tem- 
perature crossovers near quantum critical points in Section [3] In Section |U we will 
focus on the important quantum critical region and present a general discussion of 
its transport properties. An important recent development has been the complete ex- 
act solution, via the AdS/CFT correspondence, of the dynamic and transport proper- 
ties in the quantum critical region of a variety of (supersymmetric) model systems in 
two and higher dimensions: this will be described in Section|5] The exact solutions 
are found to agree with the earlier general ideas discussed here in SectionH] As has 
often been the case in the history of physics, the existence of a new class of solvable 
models leads to new and general insights which apply to a much wider class of sys- 
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tems, almost all of which are not exactly solvable. This has also been the case here, 
as we will review in Section [6] a hydrodynamic theory of the low frequency trans- 
port properties has been developed, and has led to new relations between a variety 
of thermo-electric transport co-efficients. 

The latter part of these lectures will turn to the second class of strong coupling 
problems, by describing the role of fermions near quantum critical points. In Sec- 
tion [7] we will consider some simple symmetry breaking transitions in if -wave su- 
perconductors. Such superconductors have fermionic excitations with a massless 
Dirac spectrum, and we will show how they become critical near the quantum phase 
transition. We will review how the field-theoretic 1 /N expansion does allow solu- 
tion of a large class of such problems. Finally, in Section[8]we will consider phase 
transitions of metallic systems with Fermi surfaces. We will discuss how the l/N 
expansion fails here, and review the results of recent work involving the AdS/CFT 
correspondence. 

Some portions of the discussions below have been adapted from other review 
articles by the author J2JEI. 



2 Model systems and their critical theories 
2.1 Coupled dimer antiferromagnets 

Some of the best studied examples of quantum phase transitions arise in insulators 
with unpaired 5=1/2 electronic spins residing on the sites, i, of a regular lattice. 
Using Sf (a = x,y,z) to represent the spin 5=1/2 operator on site i, the low energy 
spin excitations are described by the Heisenberg exchange Hamiltonian 



where Jy > is the antiferromagnetic exchange interaction. We will begin with a 
simple realization of this model is illustrated in Fig.[T] The 5=1/2 spins reside on 
the sites of a square lattice, and have nearest neighbor exchange equal to either J or 
J I A . Here A > 1 is a tuning parameter which induces a quantum phase transition in 
the ground state of this model. 

At A = 1, the model has full square lattice symmetry, and this case is known 
to have a Neel ground state which breaks spin rotation symmetry. This state has a 
checkerboard polarization of the spins, just as found in the classical ground state, 
and as illustrated on the left side of Fig. Q] It can be characterized by a vector order 
parameter (p a which measures the staggered spin polarization 



Hj — ^JijS" - Sj + ... 



(1) 



9 a = TJiS? 



<a 



(2) 
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Fig. 1 The coupled dimer antiferromagnet. The full red lines represent an exchange interaction J, 
while the dashed green lines have exchange J/X. The ellispes represent a singlet valence bond of 
spins (| tl) - I 4-t))A/2- 



where rj, = ±1 on the two sublattices of the square lattice. In the Neel state we have 
(<f> fl ) ^ 0, and we expect that the low energy excitations can be described by long 
wavelength fluctuations of a field (p a (x, t) over space, x, and imaginary time T. 

On the other hand, for A ^> 1 it is evident from Fig. Q] that the ground state 
preserves all symmetries of the Hamiltonian: it has total spin 5 = and can be 
considered to be a product of nearest neighbor singlet valence bonds on the J links. 
It is clear that this state cannot be smoothly connected to the Neel state, and so there 
must at least one quantum phase transition as a function A . 

Extensive quantum Monte Carlo simulations J4] |5] on this model have shown 
there is a direct phase transition between these states at a critical A c , as in Fig. Q] 
The value of A e is known accurately, as are the critical exponents characterizing 
a second-order quantum phase transition. These critical exponents are in excellent 
agreement with the simplest proposal for the critical field theory, @ which can 
be obtained via conventional Landau-Ginzburg arguments. Given the vector order 
parameter <p", we write down the action in d spatial and one time dimension, 



y w = / d d rdx 



i [{d T <p°f + v 2 {V<p°) 2 + s(<p°) 2 ] +^ [(<p«) 2 ] 2 



(3) 



as the simplest action expanded in gradients and powers of <p" which is consis- 
tent will all the symmetries of the lattice antiferromagnet. The transition is now 
tuned by varying s ~ (A — Ac). Notice that this model is identical to the Landau- 
Ginzburg theory for the thermal phase transition in a d + 1 dimensional ferromagnet, 
because time appears as just another dimension. As an example of the agreement: 
the critical exponent of the correlation length, v, has the same value, v = 0.71 1 . . ., 
to three significant digits in a quantum Monte Carlo study of the coupled dimer 
antiferromagnet, |6j and in a 5-loop analysis |7| of the renormalization group fixed 
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point of SfjjQ in d = 2. Similar excellent agreement is obtained for the double-layer 
antiferromagnet [8,9] and the coupled-plaque tte antiferromagnet. iflOl 

In experiments, the best studied realization of the coupled-dimer antiferromag- 
net is TICUCI3. In this crystal, the dimers are coupled in all three spatial dimensions, 
and the transition from the dimerized state to the Neel state can be induced by ap- 
plication of pressure. Neutron scattering experiments by Ruegg and collaborators 
ifTTI have clearly observed the transformation in the excitation spectrum across the 
transition, and these observations are in good quantitative agreement with theory] 1 ]. 



2.2 Deconfined criticality 

We now consider an analog of transition discussed in Section [27X1 but for a Hamil- 
tonian H = Hq + XH\ which has full square lattice symmetry at all X. For Hq, we 
choose a form of Hj, with /y = J for all nearest neighbor links. Thus at X = the 
ground state has Neel order, as in the left panel of Fig.Q] We now want to choose 
Hi so that increasing X leads to a spin singlet state with spin rotation symmetry re- 
stored. A large number of choices have been made in the literature, and the resulting 
ground state invariably lfl2l has valence bond solid (VBS) order; a VBS state has 
been observed in the organic antiferromagnet EtMe3P[Pd(dmit)2]2 lfl3l [141 . The 
VBS state is superficially similar to the dimer singlet state in the right panel of 
Fig. [TJ the spins primarily form valence bonds with near-neighbor sites. However, 
because of the square lattice symmetry of the Hamiltonian, a columnar arrangement 
of the valence bonds as in Fig.Q] breaks the square lattice rotation symmetry; there 
are 4 equivalent columnar states, with the valence bond columns running along dif- 
ferent directions. More generally, a VBS state is a spin singlet state, with a non-zero 
degeneracy due to a spontaneously broken lattice symmetry. Thus a direct transition 
between the Neel and VBS states involves two distinct broken symmetries: spin 
rotation symmetry, which is broken only in the Neel state, and a lattice rotation 
symmetry, which is broken only in the VBS state. The rules of Landau-Ginzburg- 
Wilson theory imply that there can be no generic second-order transition between 
such states. 

It has been argued that a second-order Neel- VBS transition can indeed occur 
fl5l . but the critical theory is not expressed directly in terms of either order parame- 
ter. It involves a fractionalized bosonic spinor Za (<* =1\ |), and an emergent gauge 
field . The key step is to express the vector field <p" in terms of Za by 

<P"=z* a KpZ P (4) 

where a" are the 2x2 Pauli matrices. Note that this mapping from (p a to Za is 
redundant. We can make a spacetime-dependent change in the phase of the Za by 
the field 0(jc,t) 

Za -► e' 6 z a (5) 
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and leave q>" unchanged. All physical properties must therefore also be invariant 
under Eq. (01, and so the quantum field theory for za has a U(l) gauge invariance, 
much like that found in quantum electrodynamics. The effective action for the za 
therefore requires introduction of an 'emergent' U(l) gauge field (where jj, =x, T 
is a three-component spacetime index). The field is unrelated the electromagnetic 
field, but is an internal field which conveniently describes the couplings between 
the spin excitations of the antiferromagnet. As we did for S^lg, we can write down 
the quantum field theory for Za and A^ by the constraints of symmetry and gauge 
invariance, which now yields 



For brevity, we have now used a "relativistically" invariant notation, and scaled 
away the spin-wave velocity v; the values of the couplings s, u are different from, but 
related to, those in S^uq. The Maxwell action for A^ is generated from short distance 
Za fluctuations, and it makes A^ a dynamical field; its coupling g is unrelated to the 
electron charge. The action 5? z is a valid description of the Neel state for s < 
(the critical upper value of s will have fluctuation corrections away from 0), where 
the gauge theory enters a Higgs phase with (z a ) 0. This description of the Neel 
state as a Higgs phase has an analogy with the Weinberg-Salam theory of the weak 
interactions — in the latter case it is hypothesized that the condensation of a Higgs 
boson gives a mass to the W and Z gauge bosons, whereas here the condensation of 
z a quenches the A^ gauge boson. As written, the s > phase of is a 'spin liquid' 
state with a S = collective gapless excitation associated with the A^ photon. Non- 
perturbative effects 0~2j associated with the monopoles in A^ (not discussed here), 
show that this spin liquid is ultimately unstable to the appearance of VBS order. 

Numerical studies of the Neel- VBS transition have focussed on a specific lattice 
antiferromagnet proposed by Sandvik lfT6l[T7l [181 . There is strong evidence for VBS 
order proximate to the Neel state, along with persuasive evidence of a second-order 
transition. However, some studies [19, 20 J support a very weak first order transition. 



2.3 Graphene 

The last few years have seen an explosion in experimental and theoretical stud- 
ies II2TI of graphene: a single hexagonal layer of carbon atoms. At the currently 
observed temperatures, there is no evident broken symmetry in the electronic exci- 
tations, and so it is not conventional to think of graphene as being in the vicinity of 
a quantum critical point. However, graphene does indeed undergo a bona fide quan- 
tum phase transition, but one without any order parameters or broken symmetry. 
This transition may be viewed as being 'topological' in character, and is associated 
with a change in nature of the Fermi surface as a function of carrier density. 




(6) 
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Pure, undoped graphene has a conical electronic dispersion spectrum at two 
points in the Brillouin zone, with the Fermi energy at the particle-hole symmet- 
ric point at the apex of the cone. So there is no Fermi surface, just a Fermi point, 
where the electronic energy vanishes, and pure graphene is a 'semi-metal'. By ap- 
plying a gate voltage, the Fermi energy can move away from this symmetric point, 
and a circular Fermi surface develops, as illustrated in Fig. [2] The Fermi surface is 




electron-like for one sign of the bias, and hole-like for the other sign. This change 
from electron to hole character as a function of gate voltage constitutes the quan- 
tum phase transition in graphene. As we will see below, with regard to its dynamic 
properties near zero bias, graphene behaves in almost all respects like a canonical 
quantum critical system. 

The field theory for graphene involves fermionic degrees of freedom. Repre- 
senting the electronic orbitals near one of the Dirac points by the two-component 
fermionic spinor W s , where s is a sublattice index (we suppress spin and 'valley' 
indices), we have the effective electronic action 

S*e = jd 2 r JdxW] [{d x + iA % -ii)8 s ^ + iv F ^d x + iv F ^dy\ %, 

where t'j are Pauli matrices in the sublattice space, pi is the chemical potential, vp 
is the Fermi velocity, and A T is the scalar potential mediating the Coulomb interac- 
tion with coupling g 1 = e 2 /e (e is a dielectric constant). This theory undergoes a 
quantum phase transition as a function of jj., at jU = 0, similar in many ways to that 
of as a function of s. The interaction between the fermionic excitations here has 
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coupling g 2 , which is the analog of the non-linearity u in S^jjg- The strength of the in- 
teractions is determined by the dimensionless 'fine structure constant' a=g 2 / {fivp) 
which is of order unity in graphene. While u flows to a non-zero fixed point value 
under the renormalization group, a flows logarithmically slowly to zero. For many 
purposes, it is safe to ignore this flow, and to set a equal to a fixed value. 



3 Finite temperature crossovers 

The previous section has described four model systems at T = 0: we examined the 
change in the nature of the ground state as a function of some tuning parameter, 
and motivated a quantum field theory which describes the low energy excitations on 
both sides of the quantum critical point. 

We now turn to the important question of the physics at non-zero temperatures. 
All of the models share some common features, which we will first explore for 
the coupled dimer antiferromagnet. For X > X c (or s > in S^uq), the excitations 
consist of a triplet of S — 1 particles (the 'triplons'), which can be understood per- 
turbatively in the large A expansion as an excited 5 = 1 state on a dimer, hopping 
between dimers (see Fig.0. The mean field theory tells us that the excitation energy 
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Fig. 3 Finite temperature crossovers of the coupled dimer antiferromagnet in Fig.fT] 



of this dimer vanishes as y/s upon approaching the quantum critical point. Fluctua- 
tions beyond mean field, described by S'i/j, show that the exponent is modified to 
s zV ', where z = 1 is the dynamic critical exponent, and v is the correlation length 
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exponent. Now imagine turning on a non-zero temperature. As long as T is smaller 
than the triplon gap, i.e. T < s zV , we expect a description in terms of a dilute gas 
of thermally excited triplon particles. This leads to the behavior shown on the right- 
hand-side of Fig. |2 delimited by the crossover indicted by the dashed line. Note that 
the crossover line approaches T = only at the quantum critical point. 

Now let us look a the complementary behavior at T > on the Neel-ordered 
side of the transition, with s < 0. In two spatial dimensions, thermal fluctuations 
prohibit the breaking of a non-Abelian symmetry at all T > 0, and so spin rotation 
symmetry is immediately restored. Nevertheless, there is an exponentially large spin 
correlation length, £, , and at distances shorter than £, we can use the ordered ground 
state to understand the nature of the excitations. Along with the spin-waves, we also 
found the longitudinal 'Higgs' mode with energy \/—2s in mean field theory. Thus, 
just as was this case for s > 0, we expect this spin-wave+Higgs picture to apply at 
all temperatures lower than the natural energy scale; i.e. for T < (— s) zV . This leads 
to the crossover boundary shown on the left-hand-side of Fig. [3] 

Having delineated the physics on the two sides of the transition, we are left 
with the crucial quantum critical region in the center of Fig. [3] This is present for 
T > |s| zV , i.e. at higher temperatures in the vicinity of the quantum critical point. To 
the left of the quantum critical region, we have a description of the dynamics and 
transport in terms of an effectively classical model of spin waves: this is the 'renor- 
malized classical' regime of Ref . |22l . To the right of the quantum critical region, we 
again have a regime of classical dynamics, but now in terms of a Boltzmann equa- 
tion for the triplon particles. A key property of quantum critical region is that there 
is no description in terms of either classical particles or classical waves at the times 
of order the typical relaxation time, T, , of thermal excitations. Instead, quantum and 
thermal effects are equally important, and involve the non-trivial dynamics of the 
fixed-point theory describing the quantum critical point. Note that while the fixed- 
point theory applies only at a single point (A = A c ) at T = 0, its influence broadens 
into the quantum critical region at T > 0. Because there is no characteristing energy 
scale associated with the fixed-point theory, kgT is the only energy scale available 
to determine X r at non-zero temperatures. Thus, in the quantum critical region ||231 

„ h 

Xr = (8) 
k B T 

where c € is a universal constant dependent only upon the universality class of the 
fixed point theory i.e. it is universal number just like the critical exponents. This 
value of x r determines the 'friction coefficients' associated with the dissipative re- 
laxation of spin fluctuations in the quantum critical region. It is also important for 
the transport co-efficients associated with conserved quantities, and this will be dis- 
cussed in Section |4] 

Let us now consider the similar T > crossovers for the other models of Sec- 
tion |2] 

The Neel-VBS transition of Section [2721 has crossovers very similar to those in 
Fig. |3] with one important difference. The VBS state breaks a discrete lattice sym- 
metry, and this symmetry remains broken for a finite range of non-zero temper- 



10 



Subir Sachdev 



atures. Thus, within the right-hand 'triplon gas' regime of Fig. [3] there is a phase 
transition line at a critical temperature 7Vbs- The value of 7Vbs vanishes very rapidly 
as s \ 0, and is controlled by the non-perturbative monopole effects which were 
briefly noted in Section l2~2l 
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Fig. 4 Finite temperature crossovers of graphene as a function of electron density n (which is tuned 
by ll in Eq. 0) and temperature, T. Adapted from Ref. 1241 . 



For graphene, the discussion above applied to Fig.|2]leads to the crossover dia- 
gram shown in Fig.|U as noted by Sheehy and Schmalian l24l . We have the Fermi 
liquid regimes of the electron- and hole-like Fermi surfaces on either side of the 
critical point, along with an intermediate quantum critical Dirac liquid. A new fea- 
ture here is related to the logarithmic flow of the dimensionless 'fine structure con- 
stant' a controlling the Coulomb interactions, which was noted in Section l2~3l In 
the quantum critical region, this constant takes the typical value a ~ l/ln(l/r). 
Consequently for the relaxation time in Eq. $Q we have ~ ln 2 (l/r). This time 
determines both the width of the electron spectral functions, and also the transport 
co-efficients, as we will see in Section |4] 



4 Quantum critical transport 

We now turn to the 'transport' properties in the quantum critical region: we con- 
sider the response functions associated with any globally conserved quantity. For 
the antiferromagnetic systems in Sections [2 . 1 1 and l2?2l this requires consideration of 
the transport of total spin, and the associated spin conductivities and diffusivities. 
For graphene, we can consider charge and momentum transport. Our discussion be- 
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low will also apply to the superfluid-insulator transition: for bosons in a periodic 
potential, this transition is described |25ll by a field theory closely related to that in 
Eq. (0). However, we will primarily use a language appropriate to charge transport 
in graphene below. We will describe the properties of a generic strongly-coupled 
quantum critical point and mention, where appropriate, the changes due to the log- 
arithmic flow of the coupling in graphene. 

In traditional condensed matter physics, transport is described by identifying 
the low-lying excitations of the quantum ground state, and writing down 'transport 
equations' for the conserved charges carried by them. Often, these excitations have 
a particle-like nature, such as the 'triplon' particles of Fig.0or the electron or hole 
quasiparticles of the Fermi liquids in Fig. [4] In other cases, the low-lying excitations 
are waves, such as the spin-waves in Fig. and their transport is described by a 
non-linear wave equation (such as the Gross-Pitaevski equation). However, as we 
have discussed in Section neither description is possible in the quantum critical 
region, because the excitations do not have a particle-like or wave-like character. 

Despite the absence of an intuitive description of the quantum critical dynam- 
ics, we can expect that the transport properties should have a universal character 
determined by the quantum field theory of the quantum critical point. In addition 
to describing single excitations, this field theory also determines the S-matrix of 
these excitations by the renormalization group fixed-point value of the couplings, 
and these should be sufficient to determine transport properties 0261 . The transport 
co-efficients, and the relaxation time to local equilibrium, are not proportional to a 
mean free scattering time between the excitations, as is the case in the Boltzmann 
theory of quasiparticles. Such a time would typically depend upon the interaction 
strength between the particles. Rather, the system behaves like a "perfect fluid" in 
which the relaxation time is as short as possible, and is determined universally by the 
absolute temperature, as indicated in Eq. (0. Indeed, it was conjectured in Ref. l27l 
that the relaxation time in Eq. (0 is a generic lower bound for interacting quan- 
tum systems. Thus the non-quantum-critical regimes of all the phase diagrams in 
Section0have relaxation times which are all longer than Eq. (0. 

The transport co-efficients of this quantum-critical perfect fluid also do not de- 
pend upon the interaction strength, and can be connected to the fundamental con- 
stants of nature. In particular, the electrical conductivity, a, is given by (in two 
spatial dimensions) |261 

a Q = e —0 a , (9) 

where <P a is a universal dimensionless constant of order unity, and we have added 
the subscript Q to emphasize that this is the conductivity for the case of graphene 
with the Fermi level at the Dirac point (for the superfluid-insulator transition, this 
would correspond to bosons at integer filling) with no impurity scattering, and at 
zero magnetic field. Here e* is the charge of the carriers: for a superfluid-insulator 
transition of Cooper pairs, we have e* = 2e, while for graphene we have e* = e. 
The renormalization group flow of the 'fine structure constant' a of graphene to 
zero at asymptotically low T, allows an exact computation in this case [28): <&a ~ 
0.05 In 2 (1 /T). For the superfluid-insulator transition, is T-independent (this is 
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the generic situation with non-zero fixed point values of the interaction |29l ) but it 
has only been computed |27l |26l to leading order in expansions in \/N (where N 
is the number of order parameter components) and in 3 — d (where d is the spatial 
dimensionality). However, both expansions are neither straightforward nor rigorous, 
and require a physically motivated resummation of the bare perturbative expansion 
to all orders. It would therefore be valuable to have exact solutions of quantum 
critical transport where the above results can be tested, and we turn to such solutions 
in the next section. 

In addition to charge transport, we can also consider momentum transport. This 
was considered in the context of applications to the quark-gluon plasma [30); appli- 
cation of the analysis of Ref. [26 1 shows that the viscosity, rj, is given by 

* = dO) 
s k B 

where s is the entropy density, and again <Pjj is a universal constant of order unity. 
The value of 0^ has recently been computed l33l for graphene, and again has a 
logarithmic T dependence because of the marginally irrelevant interaction: w 
0.008 ln 2 (l/r). 

We conclude this section by discussing some subtle aspects of the physics be- 
hind the seemingly simple result quantum-critical in Eq. For simplicity, we will 
consider the case of a "relativistically" invariant quantum critical point in 2+1 di- 
mensions (such as the field theories of Section |2~T| and |2~2l but marginally violated 
by graphene, a subtlety we ignore below). Consider the retarded correlation function 
of the charge density, %(k, (o), where k = |k| is the wavevector, and co is frequency; 
the dynamic conductivity, (j(co), is related to % by the Kubo formula, 

a{a>) = ^m^ r x{k,a>). (11) 

It was argued in Ref. (26l that despite the absence of particle-like excitations of the 
critical ground state, the central characteristic of the transport is a crossover from 
collisionless to collision-dominated transport. At high frequencies or low temper- 
atures, the limiting form for % reduces to that at T = 0, which is completely de- 
termined by relativistic and scale invariance and current conversion upto an overall 
constant 

X (k,(Q) = e —K = — , a (co) = ^-K; hco^k B T, (12) 

where K is a universal number ll3Tl . However, phase-randomizing collisions are 
intrinsically present in any strongly interacting critical point (above one spatial di- 
mension) and these lead to relaxation of perturbations to local equilibrium and the 
consequent emergence of hydrodynamic behavior. So at low frequencies, we have 
instead an Einstein relation which determines the conductivity with 
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Dk 2 „ e* 2 

X{k 1 co) = e* 2 Xc—^— , a{(o)^e* 2 x c D^—& l ® 2 ; h(o^k B T, (13) 
Dk^ — ico h 

where % c is the compressibility and D is the charge diffusion constant. Quantum 
critical scaling arguments show that the latter quantities obey 

k B T hv 2 
Xc = &» , D = Q 2 — , (14) 
k i v i kgi 

where ©\ 2 are universal numbers. A large number of papers in the literature, par- 
ticularly those on critical points in quantum Hall systems, have used the colli- 
sionless method of Eq. (fT2] i to compute the conductivity. However, the correct 
d.c. limit is given by Eq. ( fT3] l, and the universal constant in Eq. (0 is given 
by <P<j = @i@2- Given the distinct physical interpretation of the collisionless and 
collision-dominated regimes, we expect that K ^ &\©2- This has been shown in a 
resummed perturbation expansion for a number of quantum critical points |27l . 



5 Exact results for quantum critical transport 

The results of Section @] were obtained by using physical arguments to motivate 
resummations of perturbative expansions. Here we shall support the ad hoc assump- 
tions behind these results by examining an exactly solvable model of quantum crit- 
ical transport. 

The solvable model may be viewed as a generalization of the gauge theory in 
Eq. (0 to the maximal possible supersymmetry. In 2+1 dimensions, this is known 
as jY = 8 supersymmetry. Such a theory with the U(l) gauge group is free, and 
so we consider the non-Abelian Yang-Millis theory with a SU(A0 gauge group. The 
resulting supersymmetric Yang-Mills (S YM) theory has only one coupling constant, 
which is the analog of the electric charge g in Eq. (0. The matter content is naturally 
more complicated than the complex scalar Za in Eq. (O, and also involves relativistic 
Dirac fermions as in Eq. (0. However all the terms in the action for the matter fields 
are also uniquely fixed by the single coupling constant g. Under the renormalization 
group, it is believed that g flows to an attractive fixed point at a non-zero coupling 
g = g*; the fixed point then defines a supersymmetric conformal field theory in 2+1 
dimensions (a SCFT3), and we are interested here in the transport properties of this 
SCFT3. 

A remarkable recent advance has been the exact solution of this SCFT3 in the 
N — >. 00 limit using the AdS/CFT correspondence ||34| . The solution proceeds by 
a dual formulation as a four-dimensional supergravity theory on a spacetime with 
uniform negative curvature: anti-de Sitter space, or AdSzj. Remarkably, the solution 
is also easily extended to non-zero temperatures, and allows direct computation of 
the correlators of conserved charges in real time. At T > a black hole appears in 
the gravity, resulting in an AdS-Schwarzschild spacetime, and T is also the Hawking 
temperature of the black hole; the real time solutions also extend to T > 0. 
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The results of a full computation ||35l of the density correlation function, %(k, (o) 
are shown in Fig. [5] and [6] The most important feature of these results is that 
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Fig. 5 Spectral weight of the density correlation function of the SCFT3 with jY = 8 supersym- 
metry in the collisionless regime. 
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Fig. 6 As in Fig. [5] but for the collision-dominated regime. 

the expected limiting forms in the collisionless (Eq. (fT2l i) and collision-dominated 
(Eq. ( fT3] ) are obeyed. Thus the results do display the collisionless to collision- 
dominated crossover at a frequency of order kgT /%, as was postulated in Section[4] 
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An additional important feature of the solution is apparent upon describing the 
full structure of both the density and current correlations. Using spacetime indices 
(ju, V = t,x,y) we can represent these as the tensor X^v(k, co), where the previously 
considered % = % tt . At T > 0, we do not expect Xfiv to be relativistically co variant, 
and so can only constrain it by spatial isotropy and density conservation. Introduc- 
ing a spacetime momentum p^ = (o,k), and setting the velocity v = 1, these two 
constraints lead to the most general form 

^ v (k,0)) = ^-v / P I (^v^ r (^®)+^v^®)) d5) 

where p 2 = rj^p^pv with r]^ v = diag(— 1,1, 1), and P^ v and P^ v are orthogonal 
projectors defined by 

^oo = Poi = Pfa = , p Ij = ^ij^~^f > = (jlnv pr - ) - ^uvi (16) 

with the indices i,j running over the 2 spatial components. The two functions 
K T ' L (k, co) define all the correlators of the density and the current, and the results in 
Eqs. ( fl~3b and ( flZb are obtained by taking suitable limits of these functions. We will 
also need below the general identity 

K T {Q,co) =K L (0,co), (17) 

which follows from the analyticity of the T > current correlations at k = 0. 

The relations of the previous paragraph are completely general and apply to any 
theory. Specializing to the AdS-Schwarzschild solution of SYM3, the results were 
found to obey a simple and remarkable identity |35l : 

K L {k,co)K T (k,co)=J(r 2 (18) 

where J(f is a known pure number, independent of co and k. It was also shown that 
such a relation applies to any theory which is equated to classical gravity on AdSzt, 
and is a consequence of the electromagnetic self-duality of its four-dimensional 
Maxwell sector. The combination of Eqs. dTD i and (IT8b fully determines the X/iv 
correlators at k = 0: we find K L {0, co) = K T (0, co) = X, from which it follows that 
the k = conductivity is frequency independent and that <P a = &1&2 — K = J(f. 
These last features are believed to be special to theories which are equivalent to 
classical gravity, and not hold more generally. 

We can obtain further insight into the interpretation of Eq. ( fT8l ) by consider- 
ing the field theory of the superfluid-insulator transition of lattice bosons at integer 
filling. As we noted earliear, this is given by the field theory in Eq. (f3]) with the 
field <p a having 2 components. It is known that this 2-component theory of rela- 
tivistic bosons is equivalent to a dual relativistic theory, y of vortices, under the 
well-known 'particle-vortex' duality [36). Ref. ||35l considered the action of this 
particle-vortex duality on the correlation functions in Eq. < fT~5T >, and found the fol- 



16 



Subir Sachdev 



lowing interesting relations: 

K L {k,G>)K T (k,(0) = \ , K T {k,co)K L (k,co) = 1 (19) 

where K L T determine the vortex current correlations in 5? as in Eq. (fl~5T >. Unlike 
Eq. ( fTSI l. Eq. (JT9J does not fully determine the correlation functions at k = 0: it only 
serves to reduce the 4 unknown functions K L J , K L T to 2 unknown functions. The 
key property here is that while the theories S^hq and 5? are dual to each other, they 
are not equivalent, and the theory S^iq is not self-dual. 

We now see that Eq. < fT~8T > implies that the classical gravity theory of SYM3 is 
self-dual under an analog of particle-vortex duality ll35l . It is not expected that this 
self-duality will hold when quantum gravity corrections are included; equivalently, 
the SYM3 at finite N is expected to have a frequency dependence in its conductiv- 
ity at k = 0. If we apply the AdS/CFT correspondence to the superfluid-insulator 
transition, and approximate the latter theory by classical gravity on AdS4, we im- 
mediately obtain the self-dual prediction for the conductivity, <P a = 1 . This value is 
not far from that observed in numerous experiments, and we propose here that the 
AdS/CFT correspondence offers a rationale for understanding such observations. 



6 Hydrodynamic theory 

The successful comparison between the general considerations of Section@] and the 
exact solution using the AdS/CFT correspondence in Section [5] emboldens us to 
seek a more general theory of low frequency (hco <C kgT) transport in the quantum 
critical regime. We will again present our results for the special case of a relativistic 
quantum critical point in 2+1 dimensions (a CFT3), but it is clear that similar con- 
siderations apply to a wider class of systems. Thus we can envisage applications to 
the superfluid-insulator transition, and have presented scenarios under which such a 
framework can be used to interpret measurements of the Nernst effect in the cuprates 
11321 . We have also described a separate set of applications to graphene 11281 : while 
graphene is strictly not a CFT3, the Dirac spectrum of electrons leads to many simi- 
lar results, especially in the inelastic collision-dominated regime associated with the 
quantum critical region. These results on graphene are reviewed in a separate paper 
(37), where explicit microscopic computations are also discussed. 

Our idea is to relax the restricted set of conditions under which the results of Sec- 
tion|4]were obtained. We will work within the quantum critical regimes of the phase 
diagrams of Section [3] but now allow a variety of additional perturbations. First, 
we will move away from the particle-hole symmetric case, allow a finite density of 
carriers. For graphene, this means that jx is no longer pinned at zero; for the anti- 
ferromagnets, we can apply an external magnetic field; for the superfluid-insulator 
transition, the number density need not be commensurate with the underlying lat- 
tice. For charged systems, such as the superfluid-insulator transition or graphene, 
we allow application of an external magnetic field. Finally, we also allow a small 
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density of impurities which can act as a sink of the conserved total momentum of the 
CFT3. In all cases, the energy scale associated with these perturbations is assumed 
to be smaller than the dominant energy scale of the quantum critical region, which 
is kgT. The results presented below were obtained in two separate computations, 
associated with the methods described in Sections [4] and [5] and are described in the 
two subsections below. 



6.1 Relativistic magnetohydrodynamics 

With the picture of relaxation to local equilibrium at frequencies hoo <C kgT devel- 
oped in Ref. l26l . we postulate that the equations of relativistic magnetohydrody- 
namics should describe the low frequency transport. The basic principles involved in 
such a hydrodynamic computation go back to the nineteenth century: conservation 
of energy, momentum, and charge, and the constraint of the positivity of entropy 
production. Nevertheless, the required results were not obtained until our recent 
work 1321 : the general case of a CFT3 in the presence of a chemical potential, 
magnetic field, and small density of impurities is very intricate, and the guidance 
provided by the dual gravity formulation was very helpful to us. In this approach, 
we do not have quantitative knowledge of a few transport co-efficients, and this is 
complementary to our ignorance of the effective couplings in the dual gravity theory 
to be discussed in Section l6T2l 

The complete hydrodynamic analysis can be found in Ref. [32] . The analysis is 
intricate, but is mainly a straightforward adaption of the classic procedure outlined 
by Kadanoff and Martin BT1 to the relativistic field theories which describe quantum 
critical points. We list the steps: 

1. Identify the conserved quantities, which are the energy-momentum tensor, T^ v , 
and the particle number current, J 11 . 

2. Obtain the real time equations of motion, which express the conservation laws: 

dv j^ =F tiv Jv ^ d^=Q; (20) 

here F^ v represents the externally applied electric and magnetic fields which 
can change the net momentum or energy of the system, and we have not written 
a term describing momentum relaxation by impurities. 

3. Identify the state variables which characterize the local thermodynamic state — 
we choose these to be the density, p, the temperature T, and an average velocity 
it". 

4. Express T^ v and /" in terms of the state variables and their spatial and temporal 
gradients; here we use the properties of the observables under a boost by the 
velocity u^, and thermodynamic quantities like the energy density, e, and the 
pressure, P, which are determined from T and p by the equation of state of the 
CFT. We also introduce transport co-efficients associated with the gradient terms. 
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5. Express the equations of motion in terms of the state variables, and ensure that 
the entropy production rate is positive |42J. This is a key step which ensures 
relaxation to local equilibrium, and leads to important constraints on the transport 
co-efficients. In d = 2, it was found that situations with the velocity spacetime 
independent are characterized by only a single independent transport co-efficient 
11321 . This we choose to be the longitudinal conductivity at B = 0. 

6. Solve the initial value problem for the state variables using the linearized equa- 
tions of motion. 

7. Finally, translate this solution to the linear response functions, as described in 
Ref. ED. 



6.2 Dyonic black hole 

Given the success of the AdS/CFT correspondence for the specific supersymmetric 
model in Section[3] we boldly assume a similar correspondence for a generic CFT3. 
We assume that each CFT3 is dual to a strongly-coupled theory of gravity on AdSzt. 
Furthermore, given the operators associated with the perturbations away from the 
pure CFT3 we want to study, we can also deduce the corresponding perturbations 
away from the dual gravity theory. So far, this correspondence is purely formal and 
not of much practical use to us. However, we now restrict our attention to the hy- 
drodynamic, collision dominated regime, hco <C ksT of the CFT3. We would like to 
know the corresponding low energy effective theory describing the quantum gravity 
theory on AdS4. Here, we make the simplest possible assumption: the effective the- 
ory is just the Einstein-Maxwell theory of general relativity and electromagnetism 
on AdSzt. As in Section [5] the temperature T of CFT3 corresponds to introducing 
a black hole on AdS4 whose Hawking temperature is T. The chemical potential, 
H, of the CFT3 corresponds to an electric charge on the black hole, and the ap- 
plied magnetic field maps to a magnetic charge on the black hole. Such a dynoic 
black hole solution of the Einstein-Maxwell equations is, in fact, known: it is the 
Reissner-Nordstrom black hole. 

We solved the classical Einstein-Maxwell equations for linearized fluctuations 
about the metric of a dyonic black hole in a space which is asymptotically AdSzt. 
The results were used to obtain correlators of a CFT3 using the prescriptions of 
the AdS/CFT mapping. As we have noted, we have no detailed knowledge of the 
strongly-coupled quantum gravity theory which is dual to the CFT3 describing the 
superfluid-insulator transition in condensed matter systems, or of graphene. Nev- 
ertheless, given our postulate that its low energy effective field theory essentially 
captured by the Einstein-Maxwell theory, we can then obtain a powerful set of re- 
sults for CFT3s. 
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6.3 Results 



In the end, we obtained complete agreement between the two independent computa- 
tions in Sections ^. 1 l and ROl after allowing for their distinct equations of state. This 
agreement demonstrates that the assumption of a low energy Einstein-Maxwell ef- 
fective field theory for a strongly coupled theory of quantum gravity is equivalent 
to the assumption of hydrodynamic transport for Tico <C kgT in a strongly coupled 
CFT3. 

Finally, we turn to our explicit results for quantum critical transport with Tico <C 
k B T. 

First, consider adding a chemical potential, ji, to the CFT3. This will induce a 
non-zero number density of carriers p. The value of p is defined so that the total 
charge density associated with p is e*p. Then the electrical conductivity at a fre- 
quency co is 

e* 2 e* 2 pV 1 

g(m) = -^ + _ ( _ /fi)+1/Timp) (21) 

In this section, we are again using the symbol v to denote the characteristic veloc- 
ity of the CFT3 because we will need c for the physical velocity of light below. 
Here e is the energy density and P is the pressure of the CFT3. We have assumed a 
small density of impurities which lead to a momentum relaxation time Ti mp ll32l[38ll . 
In general, <P a , p, e, P, and 1 / li mp will be functions of ji/kgT which cannot be 
computed by hydrodynamic considerations alone. However, apart from <P a , these 
quantities are usually amenable to direct perturbative computations in the CFT3, or 
by quantum Monte Carlo studies. The physical interpretation of Eq. (fSTJ should be 
evident: adding a charge density p leads to an additional Drude-like contribution to 
the conductivity. This extra current cannot be relaxed by collisions between the un- 
equal density of particle and hole excitations, and so requires an impurity relaxation 
mechanism to yield a finite conductivity in the d.c. limit. 

Now consider thermal transport in a CFT3 with a non-zero fi. The d.c. thermal 
conductivity, K, is given by 

— (f)(f^) 2 . 

in the absence of impurity scattering, 1 /T; mp — » 0. This is a Wiedemann-Franz-like 
relation, connecting the thermal conductivity to the electrical conductivity in the 
jj. = CFT Note that K diverges as p — s- 0, and so the thermal conductivity of the 
jU = CFT is infinite. 

Next, turn on a small magnetic field B; we assume that B is small enough that 
the spacing between the Landau levels is not as large as kgT. The case of large 
Landau level spacing is also experimentally important, but cannot be addressed by 
the present analysis. Initially, consider the case n = 0. In this case, the result Eq. (1221 
for the thermal conductivity is replaced by 
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1 fk 2 B T\f 



<P a \ h ) \k B TB/(hc/e* 



(23) 



also in the absence of impurity scattering, 1 /Zi mp — > 0. This result relates K to the 
electrical resistance at criticality, and so can be viewed as Wiedemann-Franz-like 
relation for the vortices. A similar dependence of K appeared in the Boltzmann 
equation analysis of Ref. 091 . but our more general analysis applies in a wider and 
distinct regime 0281 . and relates the co-efficient to other observables. 

We have obtained a full set of results for the frequency-dependent thermo-electric 
response functions at non-zero B and jj,. The results are lengthy and we refer the 
reader to Ref. 021 for explicit expressions. Here we only note that the character- 
istic feature [32, 40] of these results is a new hydrodynamic cyclotron resonance. 
The usual cyclotron resonance occurs at the classical cyclotron frequency, which is 
independent of the particle density and temperature; further, in a Galilean-invariant 
system this resonance is not broadened by electron-electron interactions alone, and 
requires impurities for non-zero damping. The situtation for our hydrodynamic res- 
onance is very different. It occurs in a collision-dominated regime, and its frequency 
depends on the density and temperature: the explicit expression for the resonance 
frequency is 

e*Bpv 2 

(O c = , , DV 24 
c(e + P) 

Further, the cyclotron resonance involves particle and hole excitations moving in 
opposite directions, and collisions between them can damp the resonance even in 
the absence of impurities. Our expression for this intrinsic damping frequency is 



*2 d? 2 

e B v 

Y= -k^^Wy (25) 

relating it to the quantum-critical conductivity as a measure of collisions between 
counter-propagating particles and holes. We refer the reader to a separate discus- 
sion |28l of the experimental conditions under which this hydrodynamic cyclotron 
resonance may be observed. 



7 d- wave superconductors 

We now turn to the second class of strong-coupling problems outlined in Section[TJ 
those involving quantum critical points with fermionic excitations. This section will 
consider the simpler class of problems in which the fermions have a Dirac spectrum, 
and the field-theoretic 1 /N expansion does allow for substantial progress. 

We will begin in Section l7TT1 bv an elementary discussion of the origin of these 
Dirac fermions. Then we will consider two quantum phase transitions, both involv- 
ing a simple Ising order parameter. The first in Section lT2l with time-reversal sym- 
metry breaking, leads to a relativistic quantum field theory closely related to the 
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Gross-Neveu model. The second model of Section [731 involves breaking of a lattice 
rotation symmetry, leading to "Ising-nematic" order. The theory for this model is not 
relativistically invariant: it is strongly coupled, but can be controlled by a traditional 
\/N expansion. 

We note that symmetry breaking transitions in graphene are also described by 
field theories similar to those discussed in this section 143 . 1441 . 



7.1 Dirac fermions 

We begin with a review of the standard BCS mean-field theory for a cZ-wave su- 
perconductor on the square lattice, with an eye towards identifying the fermionic 
Bogoliubov quasiparticle excitations. For now, we assume we are far from any QPT 
associated with SDW, Ising-nematic, or other broken symmetries. We consider the 
Hamiltonian 

Hu = £ e k c\ a c ka +7i £ Si ■ Sj (26) 

* m 

where cj a is the annihilation operator for an electron on site j with spin a =t, I, c ka 
is its Fourier transform to momentum space, £ k is the dispersion of the electrons (it 
is conventional to choose E k = — 2t\ (cos(k x ) + cos(fcy)) — 2t2{cos{k x + k y ) + cos{k x — 
k y )) — /i, with t\ 2 the first/second neighbor hopping and jj. the chemical potential), 
and the J\ term is similar to that in Eq. (HJ with 



St, 



(27) 



and a" the Pauli matrices. We will consider the consequences of the further neighbor 
exchange interactions for the superconductor in Section 17.21 below. Applying the 
BCS mean-field decoupling to H t j we obtain the Bogoliubov Hamiltonian 



Hbcs =Y< e k c la c k<* - T E 4 M 



JH 



h.c. 



(28) 



For a wide range of parameters, the ground state energy optimized by a d x 2_ y 2 wave- 
function for the Cooper pairs: this corresponds to the choice A x 



-A y = A x 2_ y i. 



The value of A x 2_ y 2 is determined by minimizing the energy of the BCS state 



d 2 k 
An 2 



(29) 



where the fermionic quasiparticle dispersion is 



E t = 



J\A x 2_ y 2 {cosk x — cosky) 



1/2 



(30) 
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The energy of the quasiparticles, E^, vanishes at the four points (±Q,±Q) at 
which Ek = 0. We are especially interested in the low energy quasiparticles in the 
vicinity of these points, and so we perform a gradient expansion of Hgcs near eacn 
of them. We label the points Qi = (Q,Q), Qi = (~Q,Q), & = (~Q,-Q), Qa = 
(Q, —Q) and write 



-ja 



■ fla(n)e^ +f2a(?j)e&- r >+f 3 aVj)&- fl +f4a(? J )e&'j, 



(31) 



while assuming the /i-4,a(?) are slowly varying functions of x. We also introduce 
the bispinors fi = (/i t ,/ 3 \,/u, -/ 3 f t ), and f 2 = (Z^,/^,/^-/^). and then ex- 
press Hgcs m terms of f j 2 while performing a spatial gradient expansion. This 
yields the following effective action for the fermionic quasiparticles: 



dxd l r 



W2 



V2 K 



(32) 



where the T x ' z are 4x4 matrices which are block diagonal, the blocks consisting of 
2x2 Pauli matrices. The velocities vf.a are given by the conical structure of near 
the Q1-4. we have vp = \^k£k\k=Q a I and = |7i4 r 2_ Y 2V2sin(g)|. In this limit, the 
energy of the W\ fermionic excitations is E\ = (vp (k x + k y ) 2 /2 + v\ (k x — k y ) 2 /2) l l 2 
(and similarly for f^X which is the spectrum of massless Dirac fermions. 



7.2 Time-reversal symmetry breaking 

We will consider a simple model in which the pairing symmetry of the supercon- 
ductor changes from d x 2_ y i to d x i_ y i ± idxy. The choice of the phase between the 
two pairing components leads to a breaking of time-reversal symmetry. Studies of 
this transition were originally motivated by the cuprate phenomenology, but we will 
not explore this experimental connection here because the evidence has remained 
sparse. 

The mean field theory of this transition can be explored entirely within the con- 
text of BCS theory, as we will review below. However, fluctuations about the BCS 
theory are strong, and lead to non-trivial critical behavior involving both the col- 
lective order parameter and the Bogoliubov fermions: this is probably the earliest 
known example 1431 l46l of the failure of BCS theory in two (or higher) dimensions 
in a superconducting ground state. At T > 0, this failure broadens into the "quantum 
critical" region. 

We extend H t j in Eq. d26l i so that BCS mean-field theory permits a region with 
dxy superconductivity. With a J2 second neighbor interaction, Eq. d26l i is modified 
to: 
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We will follow the evolution of the ground state of H t j as a function of J%lJ\. 
The mean-field Hamiltonian is now modified from Eq. (l28l l to 

Hues = L £ kc\a c ko - y E 4 m ( c ]t c J+A,4 - C ]l c ]+A,t) + h ' c ' 
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(33) 



h 



h.c. 



(34) 



where the second summation over v is along the diagonal neighbors x+y and — i + 
y. To obtain of XY pairing along the diagonals, we choose A x +y = —A- x +y = A_ xy . We 
summarize our choices for the spatial structure of the pairing amplitudes (which 
determine the Cooper pair wavefunction) in Fig [7] The values of A x 2_ y 2 and Axy are 




x2-y2 



A x 2. y2 A xy 



Fig. 7 Values of the pairing amplitudes, —(c^Cji — c,i Cyf) with i the central site, and j is one of 
its 8 near neighbors. 



to be determined by minimizing the ground state energy (generalizing Eq. < f29b > 

d 2 k 



where the quasiparticle dispersion is now (generalizing Eq. ( |30l >) 

J\ A x 2 2 (cos k x — cos fey) + 2/2 sin £ v sin k y 



-HA^ 



[Ek - e k ] 



1/2 



(35) 



(36) 



Notice that the energy depends upon the relative phase of A x 2_ y 2 and A^r. this phase 
is therefore an observable property of the ground state. 

It is a simple matter to numerically carry out the minimization of Eq. 1361 . and 
the results for a typical choice of parameters are shown in Fig[8]as a function J%jJ\. 
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Fig. 8 BCS solution of the phenomenological Hamiltonian H,j in Eq. j33> . Shown are the opti- 
mum values of the pairing amplitudes % | and \A xy \ as a function of J% for t\ = 1, ?2 = —0.25, 
fi = —1.25, and J\ fixed at J\ = 0.4. The relative phase of the pairing amplitudes was always 
found to obey Eq. d37t . The dashed lines denote locations of phase transitions between d x i„i, 
d x i_yi + id„, and d„ superconductors. The pairing amplitudes vanishes linearly at the first transi- 
tion corresponding to the exponent figcs = 1 in Eq. <40K The Brillouin zone location of the gapless 
Dirac points in the dg, 2 superconductor is indicated by filled circles. For the dispersion e^. ap- 
propriate to the cuprates, the d xy superconductor is fully gapped, and so the second transition is 
ordinary Ising. 



One of the two amplitudes A x 2_ y 2 or A xy is always non-zero and so the ground 
state is always superconducting. The transition from pure d x 2_ y 2 superconductivity 
to pure d xy superconductivity occurs via an intermediate phase in which both order 
parameters are non-zero. Furthermore, in this regime, their relative phase is found 
to be pinned to ±7r/2 i.e. 



arg(Arv) = arg(A x 2_ y 2)± n/2 



(37) 



The reason for this pinning can be intuitively seen from Eq. d36t : only for these 
values of the relative phase does the equation = never have a solution. In other 
words, the gapless nodal quasiparticles of the d x 2_ y 2 superconductor acquire a finite 
energy gap when a secondary pairing with relative phase ±n/2 develops. By a level 
repulsion picture, we can expect that gapping out the low energy excitations should 
help lower the energy of the ground state. The intermediate phase obeying Eq. d37l i 
is called a d x 2_ y 2 + idxy superconductor. 

The choice of the sign in Eq. d37l i leads to an overall two-fold degeneracy in the 
choice of the wavefunction for the d x 2_ y 2 + idxy superconductor. This choice is re- 
lated to the breaking of time-reversal symmetry, and implies that the d x 2_ y 2 + idxy 
phase is characterized by the non-zero expectation value of a Z2 Ising order param- 
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eter; the expectation value of this order vanishes in the two phases (the d x 2_ v 2 and 
dxy superconductors) on either side of the d x 2_ v i + id xy superconductor. As is con- 
ventional, we will represent the Ising order by a real scalar field 0. Fluctuations of 
(j> become critical near both of the phase boundaries in Fig [8] As we will explain 
below, the critical theory of the d x i_ji to d x i_ y i + idxy transition is not the usual 
4 field theory which describes the ordinary Ising transition in three spacetime di- 
mensions. (For the dispersion appropriate to the cuprates, the d xy superconductor 
is fully gapped, and so the d x i_ y i + idxy to d XJ transition in Fig. [8] will be ordinary 
Ising.) 

Near the phase boundary from d x 2_ y 2 to d x i_ y i + id xy superconductivity it is clear 
that we can identify 

$ = iA X y, (38) 

(in the gauge where A x 2_ji is real). We can now expand Ebcs m Eq. d35l l for small 
(with A x 2_ y 2 finite) and find a series with the structure P71I481 

E BCS = Eo + s<p 2 + v\^ + ..., (39) 

where s, v are coefficients and the ellipses represent regular higher order terms in 
even powers of 0; s can have either sign, whereas v is always positive. Notice the 
non-analytic \ <j) | 3 term that appears in the BCS theory — this arises from an infrared 
singularity in the integral in Eq. d35l ) over at the four nodal points of the d x 2_J2 
superconductor, and is a preliminary indication that the transition differs from that 
in the ordinary Ising model, and that the Dirac fermions play a central role. We can 
optimize by minimizing Ebcs m Eq. ( 1391 — this shows that (0) = for s > 0, and 
(0) ^ for s < 0. So s ~ (Jz/Ji)c — Ji/J\ where {J2/J\) c is the first critical value 
in Fig [8] Near this critical point, we find 

W)~[s c -sf, (40) 

where we have allowed for the fact that fluctuation corrections will shift the critical 
point from s = to s = s c . The present BCS theory yields the exponent fines = 1; 
this differs from the usual mean-field exponent Pmf = 1 /2, and this is of course due 
to the non-analytic |0| 3 term in Eq. d39l . 

We can now write down the required field theory of the onset of d xy order. In 
addition to the order parameter 0, the field theory should also involve the low energy 
nodal fermions of the d x 2_ y 2 superconductor, as described by ,5V in Eq. d32t . For 
the <p fluctuations, we write down the usual terms permitted near a phase transition 
with Ising symmetry: 



d 2 rdz 



1 -((d z $) 2 + c 2 (d x $) 2 +c 2 (d y 4) 2 +s$ 2 ) + ^ 4 



(41) 



Note that, unlike Eq. (|39l , we do not have any non-analytic |0| 3 terms in the action: 
this is because we have not integrated out the low energy Dirac fermions, and the 
terms in Eq. (f4TT > are viewed as arising from high energy fermions away from the 
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nodal points. Finally, we need to couple the <j) and f i 2 excitations. Their coupling 
is already contained in the last term in Eq. (1341 : expressing this in terms of the fi 2 
fermions using Eq. (f3TT > we obtain 



= / d z rdT 



(42) 



AfjVfi - W^x y W^ 
where / d- X y is a coupling constant. The partition function of the full theory is now 

% M = [ ®$ Wi W 2 exp (-yip -y^- , (43) 



where y^ was in Eq. ( 1321 ). It can now be checked that if we integrate out the W\ 2 
fermions for a spacetime independent <j), we do indeed obtain a |0| 3 term in the 
effective potential for (j>. 

We begin our analysis of by assuming that the transition is described by 
a fixed point with #, T = 0: then the theory for the transition would be the ordi- 
nary 4 field theory y&, and the nodal fermions would be innocent spectators. The 
scaling dimension of at such a fixed point is (1 + rj/) /2 (where rj/ is the anoma- 
lous order parameter exponent at the critical point of the ordinary three dimensional 
Ising model), while that of fi .2 is 1. Consequently, the scaling dimension of i?^, is 
(1 — T]i)/2 > 0. This positive scaling dimension implies that is relevant and the 
•&xy = fixed point is unstable: the Dirac fermions are fully involved in the critical 
theory. 

Determining the correct critical behavior now requires a full renormalization 
group analysis of Z^ui- This has been described in some detail in Ref. l46l , and 
we will not reproduce the details here. The main result we need for our purposes 
is that couplings d- xy , u, Vf/c and v^/c all reach non-zero fixed point values which 
define a critical point in a new universality class. These fixed point values, and 
the corresponding critical exponents, can be determined in expansions in either 
(3 — d) 1 45 , 46 1 (where d is the spatial dimensionality) or 1 /N |49) (where N is the 
number of fermion species). An important simplifying feature here is that the fixed 
point is actually relativistically invariant. Indeed the fixed point has the structure of 
the so-called Higgs-Yukawa (or Gross-Neveu) model which has been studied exten- 
sively in the particle physics literature f50l in a different physical context: quantum 
Monte Carlo simulation of this model also exist IfSTI , and provide probably the most 
accurate estimate of the exponents. 

The non-trivial fixed point has strong implications for the correlations of the 
Bogoliubov fermions. The fermion correlation function G\ = (fif/ ) obeys 

n n \ co + v F k x i: z + v A T x 

l( ' a>= (^+^-^)('-# ^ 

at low frequencies for s > s c . Away from the critical point in the d x 2_ y 2 supercon- 
ductor with s > s c , Eq d44b holds with r/y = 0, and this is the BCS result, with sharp 
quasi-particle poles in the Green's function. At the critical point s = s c Eq. ( PPfl i 
holds with the fixed point values for the velocities (which satisfy \'f = V4 = c) and 



Condensed matter and AdS/CFT 



27 



with the anomalous dimension r\f 7^ — the (3 — d) expansion P31 estimate is 
r\f ~ (3 — d)/l4, and the l/N expansion estimate |49j is r\f 1/(37T 2 A^), with 
N = 2. This is clearly non-BCS behavior, and the fermionic quasiparticle pole in the 
spectral function has been replaced by a branch-cut representing the continuum of 
critical excitations. The corrections to BCS extend also to correlations of the Ising 
order 0: its expectation value vanishes as Eq. d40b with the Monte Carlo estimate 
/3 ~ 0.877 IBTl . The critical point correlators of <p have the anomalous dimension 
rj « 0.754 BTI . which is clearly different from the very small value of the exponent 
rj] at the unstable = fixed point. The value of j3 is related to 77 by the usual 
scaling law j3 = (1 + rj)v/2, with v « 1.00 the correlation length exponent (which 
also differs from the exponent V/ of the Ising model). 



7.3 Nematic ordering 

We now consider an Ising transition associated with "Ising-nematic" ordering in 
the of-wave superconductor. This is associated with a spontaneous reduction of the 
lattice symmetry of the Hamiltonian from "square" to "rectangular". Our study is 
motivated by experimental observations of such a symmetry breaking in the cuprate 
superconductors ll53l l54l l52l . 

The ingredients of such an ordering are actually already present in our simple 
review of BCS theory in Section \T?\\ In Eq. (l28l ). we introduce 2 variational pair- 
ing amplitudes A x and A y . Subsequently, we assumed that the minimization of the 
energy led to a solution with d x i_ v i pairing symmetry with A x = —Ay = A x 2_ v 2- 
However, it is possible that upon including the full details of the microscopic in- 
teractions we are led to a minimum where the optimal solution also has a small 
amount of s-wave pairing. Then \A X \ 7^ |A,,|, and we can expect all physical proper- 
ties to have distinct dependencies on the x and y co-ordinates. Thus, one measure of 
the the Ising nematic order parameter is |At| 2 — |Ay| 2 . 

The derivation of the field theory for this transition follows closely our presen- 
tation in Section 17.21 We allow for small Ising-nematic ordering by introducing a 
scalar field (j) and writing 

A x = A x 2_ y 2 + (p ; A y = -A x z_ y2 + <l>. (45) 

The evolutions of the Dirac fermion spectrum under such a change is indicated in 
Fig. [9] We now develop an effective action for <p and the Dirac fermions fi 2- The 
result is essentially identical to that in Section l772l apart from a change in the struc- 
ture of the Yukawa coupling. Thus we obtain a theory S%i + ^ + ^if^ , defined by 
Eqs. d32b and (f4~Tb . and where Eq. d42l is now replaced by 



fv* = $ij d 2 rdx <j) (VjVf! + f 2 Vf> 2 



(46) 
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Fig. 9 Phase diagram of Ising nematic ordering in a rf-wave superconductor as a function of the 
coupling s in ,9'§ . The filled circles indicate the location of the gapless fermionic excitations in the 
Brillouin zone. The two choices for s < s c are selected by the sign of {(j>). 



The seemingly innocuous change between Eqs. d42l i and d46b however has strong 
consequences. This is partly linked to the fact with S^^^ cannot be relativistically 
invariant even after all velocities are adjusted to equal. A weak-coupling renormal- 
ization group analysis in powers of the coupling #/ was performed in (3 — d) di- 
mensions in Refs. B31I461 . and led to flows to strong coupling with no accessible 
fixed point: thus no firm conclusions on the nature of the critical theory were drawn. 

This problem remained unsolved until the recent works of Refs. 1 55 , 56) . It is es- 
sential that there not be any expansion in powers of the coupling This is because 
it leads to strongly non-analytic changes in the structure of the propagator, which 
have to be included at all stages. In a model with fermion flavors, the 1 /N expan- 
sion does avoid any expansion in The renormalization group analysis has to be 
carried out within the context of the 1 /N expansion, and this involves some rather 
technical analysis which is explained in Ref. (56 1. In the end, an asymptotically ex- 
act description of the vicinity of the critical point was obtained. It was found that the 
velocity ratio vf / diverged logarithmically with energy scale, leading to strongly 
anisotropic 'arc-like' spectra for the Dirac fermions. Associated singularities in the 
thermal conductivity have also been computed f57l . 



8 Metals 



This section considers symmetry breaking transitions in two-dimensional metals. 
Away from the quantum critical point, the phases will be ordinary Fermi liquids. We 
will be interested in the manner in which the Fermi liquid behavior breaks down at 
the quantum critical point. Our focus will be exclusively on two spatial dimensions: 
quantum phase transitions of metals in three dimensions are usually simpler, and the 
traditional perturbative theory appears under control. 
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In Section|7]the fermionic excitations had vanishing energy only at isolated nodal 
points in the Brillouin zone: see Figs. [8] and [9] Metals have fermionic excitations 
with vanishing energy along a line in the Brillouin zone. Thus we can expect them to 
have an even stronger effect on the critical theory. This will indeed be the case, and 
we will be led to problems with a far more complex structure. Unlike the situation in 
insulators and <i-wave superconductor, many basic issues associated with ordering 
transition in two dimensional metals have not been fully resolved. The problem 
remains one of active research and is being addressed by many different approaches. 
In recent papers l58ll59l , Metlitski and the author have argued that the problem is 
strongly coupled, and proposed field theories and scaling structures for the vicinity 
of the critical point. We will review the main ingredients for the transition involving 
Ising-nematic ordering in a metal. Thus the symmetry breaking will be just as in 
Section r7~3l but the fermionic spectrum will be quite different. Our study here is also 
motivated by experimental observations in the cuprate superconductors Il53ll54ll52ll . 

As in Section|7] let us begin by a description of the non-critical fermionic sector, 
before its coupling to the order parameter fluctuations. We use the band structure 
describing the cuprates in the over-doped region, well away from the Mott insulator. 
Here the electrons cj_ a are described by the kinetic energy in Eq. d26l i. which we 
write in the following action 



As in Section 17.31 we will have an Ising order parameter represented by the real 
scalar field 0, which is described as before by in Eq. (l4~Tb . Its coupling to the 
electrons can be deduced by symmetry considerations, and the most natural coupling 
(the analog of Eqs. d46l l) is 



where V is the volume. The momentum dependent form factor is the simplest choice 
which changes sign under x -s-> y, as is required by the symmetry properties of (j). The 
sum over q is over small momenta, while that over k extends over the entire Brillouin 
zone. The theory for the nematic ordering transition is now described by S^ c + S^q + 
5f cl h . The phase diagram as a function of the coupling s in and temperature T is 
shown in Fig.[lO] Note that there is a line of Ising phase transitions at T — T c : this 
transition is in the same universality class as the classical two-dimensional Ising 
model. However, quantum effects and fermionic excitations are crucial at T = 
critical point at s = s c and its associated quantum critical region. 

A key property of Eq. d48b is that small momentum critical fluctuations can 
efficiently scatter fermions at every point on the Fermi surface. Thus the non-Fermi 
singularities in the fermion Green's function will extend to all points on the Fermi 
surface. This behavior is dramatically different from all the field theories we have 
met so far, all of which had singularities only at isolated points in momentum space. 




(47) 




(48) 
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Fig. 10 Phase diagram of Ising nematic ordering in a metal as a function of the coupling s in 
S^a and temperature T. The Fermi surface for s > is as in the overdoped region of the cuprates, 
with the shaded region indicating the occupied hole (or empty electron) states. The choice between 
the two quadrapolar distortions of the Fermi surface is determined by the sign of (<j>). The line of 
T > phase transitions at T c is described by Onsager's solution of the classical two-dimensional 
Ising model. We are interested here in the quantum critical point at s = s c , which controls the 
quantum-critical region. 



We evidently have to write down a long-wavelength theory which has singularities 
along a line in momentum space. 

We describe the construction of Ref. l58l of a field theory with this unusual 
property. Pick a fluctuation of the order parameter <p at a momentum q. As shown 
in Fig. QT|this fluctuation will couple most efficiently to fermions near two points 
on the Fermi surface, where the tangent to the Fermi surface is parallel to q. A 
fermion absorbing momentum q at these points, changes its energy only by ~ q ; 
at all other points on the Fermi surface the change is ~ q. Thus we are led to focus 
on different points on the Fermi surface for each direction of q. In the continuum 
limit, we will therefore need a separate field theory for each pair of points ±ko on 
the Fermi surface. We conclude that the quantum critical point is described by an 
infinite number of field theories. 

There have been earlier descriptions of Fermi surfaces by an infinite number of 
field theories [ 60 6 1 62 63 64 65 66 1 . However many of these earlier works dif- 
fer in a crucial respect from the the theory to be presented here. They focus on the 
motion of fermions transverse to Fermi surface, and so represent each Fermi sur- 
face point by a 1+1 dimensional chiral fermion. Thus they have an infinite number 
of 1+1 dimensional field theories, labeled by points on the one dimensional Fermi 
surface. The original problem was 2+1 dimensional, and so this conserves the total 
dimensionality and the number of degrees of freedom. However, we have already 
argued above that the dominant fluctuations of the fermions are in a direction trans- 
verse to the Fermi surface, and so we believe this earlier approach is not suited for 




the vicinity of the nematic quantum critical point. As we will see in Section [8Tl we 
will need an infinite number of 2+ 1 dimensional field theories labeled by points on 
the Fermi surface. (Some of the earlier works included fluctuations transverse to the 
Fermi surface [65, 66), but did not account for the curvature of the Fermi surface; 
it is important to take the scaling limit at fixed curvature, as we will see.) Thus we 
have an emergent dimension and a redundant description of the degrees of freedom. 
We will see in Section liOl how compatibility conditions ensure that the redundancy 
does not lead to any inconsistencies. The emergent dimensionality suggests a con- 
nection to the AdS/CFT correspondence, as will be discussed in Section [831 



8.1 Field theories 

Let us now focus on the vicinity of the points ±&o, by introducing fermionic field 
V± by 

V+®=%+* , W-{k)=c_ lo+I . (49) 

Then we expand all terms in 5? c + J^a + S^ c a in spatial and temporal gradients. 
Using the co-ordinate system illustrated in Fig.QT| performing appropriate rescaling 
of co-ordinates, and dropping terms which can later be easily shown to be irrelevant, 
we obtain the 2+ 1 dimensional Lagrangian 

% = Via (C<5t - id x - dfj y+a + W-a (C^t + id x - d^j 1/A_ a 
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Here £, X and s are coupling constants, with s the tuning parameter across the tran- 
sition; we will see that all couplings apart from s can be scaled away or set equal 
to unity. We now allow the spin index a = 1 . . .N, as we will be interested in the 
structure of the large N expansion. Note that Eq. ( l50b has the same basic structure as 
the models considered in Section [7] apart from differences in the spatial gradients 
and the matrix structure. We will see that these seemingly minor differences will 
completely change the physical properties and the nature of the large N expansion. 



8.2 Symmetries 

A first crucial property of Jz? is that the fermion Green's functions do indeed have 
singularities along a line in momentum space, as was required by our discussion 
above. This singularity is a consequence of the invariance of Jz? under the following 
transformation 

<j)(x,y) -> <j)(x,y + dx), \j/ s (x,y) -> e- is ^+^ x) xi/ s {x,y + Ox), (51) 

where 9 is a constant. Here we have dropped time and a indices because they play 
no role, and s = ±. We can view this transformation as one which performs a 'rota- 
tion' of spatial co-ordinates, moving the point Icq to neighboring points on the Fermi 
surface. We are not assuming the Fermi surface is circular, and so the underlying 
model is not rotationally invariant. However, we are considering a limiting case 
of a rotation, precisely analogous to the manner in which Galilean transformations 
emerge as a limiting case of a relativistic transformation (x behaves like time, and y 
as space, in this analogy). This 'Galilean' symmetry is an emergent symmetry of Jz? 
for arbitrary shapes of the Fermi surface. It is not difficult to now show from dBll l 
that the Green's function D, and the fermion Green's functions G s obey the exact 
identities 

D(q x ,q y ) = D(q y ) (52) 
G s (q.x , q y ) = G(sq x + q* ) . (53) 

So we see that the f + Green's function depends only on q x + q^,. The singularities 
of this function appear when q x + q^ = 0, and this is nothing but the equation of 

the Fermi surface passing through the point Icq in Figs. QT|and[T2] Thus we have 
established the existence of a line of singularities in momentum space. 

The identities in Eq. (IBlT i also help establish the consistency of our description in 
terms of an inifinite number of 2+ 1 dimensional field theories. Consider the fermion 
Green's function at the point P in Fig. Q~2] This can be computed in terms of the 
2+1 dimensional field theory defined at the point Icq, or from that at a neighboring 
point k\. Eq. d53l l ensures that both methods yield the same result. A little geometry 
[58 1 shows that q x + q\ is an invariant that measures the distance between P and the 
closest point on the Fermi surface: thus it takes the same value in the co-ordinates 
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1 h; — yx 

Fig. 12 The fermion correlator at the point P can be described either in terms of the 2+1 dimen- 
sional field theory at X'o, or that at k\ . 

systems at Icq and k\, with q x + cfi = q' x + g' 2 , establishing the identity of the two 
computations. 



8.3 Scaling theory 

Let us now discuss the behavior of _Sf under renormalization group scaling transfor- 
mations. The structure of the spatial gradient terms in the Lagrangian indicates that 
the rescaling of spatial co-ordinates should be defined by 

x'=x/b 2 , y'=y/b. (54) 

The invariance in Eq. d53l l implies that these scalings are exact, and the spatial 
anisotropy acquires no fluctuation corrections. Or, in other words 

dim|j]=-l , dim[jc] = -2 (55) 

For now, let us keep the rescaling of the temporal co-ordinate general: 

dimfr] = -z. (56) 

Note that the dynamic critical exponent z is defined relative to the spatial co-ordinate 
y tangent to the Fermi surface (other investigators sometimes define it relative to the 
co-ordinate x normal to the Fermi surface, leading to a difference in a factor of 
2). We define the engineering dimensions of the fields so that co-efficients of the y 
derivatives remain constant. Allowing for anomalous dimensions tj^ and T]y from 
loop effects we have 



dim[0] = (l +2 + 77^/2 , dimM = (1 +z + rj v )/2. 



(57) 
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Using these transformations, we can examine the scaling dimensions of the cou- 
plings in Jz? at tree level 

dim[C]=2-z-77 v , dim[A] = (3-z-7?^-2Tj v )/2. (58) 

We will see in Section 18.41 that low order loop computations suggest that the 
anomalous dimensions rj^ and T]y are small, and that z«3. Assuming these esti- 
mates are approximately correct, we see that the coupling £ is strongly irrelevant. 
Thus we can send £ — > in all our computations. However, we do not set £ = 
at the outset, because the temporal derivative term is needed to define the proper 
analytic structure of the frequency loop integrals (67). As we will discuss later, the 
limit £ —> 0, also dramatically changes the counting of powers of 1 /N in the loop 
expansion l67ll . 

Also note that these estimates of the scaling dimensions imply dim[A] rj 0. Thus 
the fermion and order parameter fluctuations remain strongly coupled at all scales. 
Conversely, we can also say that the requirement of working in a theory with fixed 
A implies that z~3; this circumvents the appeal to loop computations for taking 
the £ — > limit. With a near zero scaling dimension for A, we cannot expand per- 
turbatively in powers of A. This features was also found in Section l731 but there we 
were able to use the \/N expansion to circumvent this problem. 

Moving beyond tree level considerations, we note that another Ward identity 
obeyed by the theory Jz? allows us to fix the scaling dimension of exactly. This 
Ward identity is linked to the fact that <p appears in the Yukawa coupling like the 
x component of a gauge field coupled to the fermions l58l . The usual arguments 
associated with gauge invariance then imply that dim[0] = 2 (the same as the scaling 
dimension of d x ), and that we can work in theory in which the "gauge coupling" A 
set equal to unity at all scales. Note that with this scaling dimension, we have the 
exact relation 

ty=3-z. (59) 

Note also that Eq. $5% now implies that dim [A] = rj v at tree level, which is the same 
as the tree level transformation of the spatial derivative terms. The latter terms have 
been set equal to unity by rescaling the fermion field, and so it is also consistent to 
set A = 1 from now on. 

We reach the remarkable conclusion that at the critical point s = s c , Jz? is in- 
dependent of all coupling constants. The only parameter left is N, and we have no 
choice but to expand correlators in powers of 1 /N. The characterization of the crit- 
ical behavior only requires computations of the exponents z and 77^, and associated 
scaling functions. 

We can combine all the above results into scaling forms for the and f Green's 
functions at the quantum critical point at T = 0. These are, respectively 
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D-\q Xi q y ,m) = qf l ^ D (^j (60) 

G-\q x ,q y ,(0) = {q x + ql) l -^l 2 & G { . "..J , (61) 



where &d and &q are scaling functions. 



8.4 Large N expansion 

We have come as far as possible by symmetry and scaling analyses alone on ££ . 
Further results require specific computations of loop corrections, and these can only 
be carried out within the context of the l/N expansion. At N = °°, the propagator 
at criticality is 

for imaginary frequencies CO. This is clearly compatible with Eq. < f60b with z = 
3. The leading correction to the fermion propagator comes from the self energy 
associated with one <p exchange, and this leads to 

G; 1 =-/Cco + ? , + ^-/sgn(co)- 7 =^^|co| 2 / 3 , (63) 

which is also compatible with Eq. (|6H with z = 3 and r\ v = 0. Notice also that as 
CO — 5-0, the £co term in Eq. (l63l is smaller than the |o| 2//3 term arising from the 
self energy at order 1 /N; this relationship is equivalent to our earlier claim that £ is 
irrelevant at long scales, and so we should take the limit £ — > to obtain our leading 
critical scaling functions. 

The structure of Eq. d63l also illustrates a key difficulty associated with the £ — >• 
limit. At £ = 0, the leading co dependence of G _1 is ~ |o| 2 / 3 /A^. Feynman graphs 
which are sensitive to this co dependence will therefore acquire additional factors 
of N, leading to a breakdown of the conventional counting of powers of 1 /N in the 
higher loop graphs. 

This breakdown of the l/N expansion was investigated by Lee ll67l for a 'sin- 
gle patch' theory with fermions only at ko (and not at —ko). We see from Eq. d63l 
that the 1 /N term in G~' becomes important when q x = q y = i.e. the fermion is 
precisely on the Fermi surface. Thus the power of N is maximized when fermions 
in all internal lines are on the Fermi surface. Such a Fermi surface restriction is sat- 
isfied only in a subspace of reduced dimension in the momentum space integral of 
any Feynman graph. Lee presented an algorithm for computing the dimensionality 
of this restricted subspace: he demonstrated that the power of 1 /N was determined 
by the genus of the surface obtained after drawing the graph in a double-line rep- 
resentation. So determining the leading l/N terms in Eqs. d62l and d63l requires 
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summation of the infinite set of planar graphs. This problem remains unsolved, but 
the unexpected appearance of planar graphs does suggest an important role for the 
AdS/CFT correspondence. 

The structure of the loop expansion for the 'two patch' theory with fermions at 
±Icq, as written in Eq. ( f50l l. was studied in Ref. l58l . It was found that z — 3 was 
preserved upto three loops, but a small non-zero value for rjy did appear at three 
loop order. Also, the genus counting of powers of 1 /N was found to break down, 
with larger powers of N appearing in some three loop graphs. 



8.5 AdS/CFT correspondence 

There has been a great deal of recent work ll^l^l70ll7Tll72ll73ll7^|75ll76ll77l 
EH [79] [80] [81] E2 EH ED investigating the structure of Fermi surfaces using the 
AdS/CFT correspondence. The results obtained so far do have features that resem- 
ble our results above for the Ising-nematic transition in two dimensional metals. 
However, the precise connection remains obscure, and is an important topic for fu- 
ture research. In particular, a microscopic understanding of the field content of the 
CFT dual of the AdS theory is lacking, although there has been interesting progress 
very recently ||79ll83l . 

One of the main results of the analysis of Ref. f70l is the general structure of the 
fermion Green's function obtained in a theory dual to a Reissner-Nordstrom black 
hole in AdS4. This had the form 

G~ l (k,(o) = -iCQ + v F (\k\~k F )-cia> 6 , (64) 

where the momentum k is now measured from the origin of momentum space (and 
not from a Fermi surface), and the complex number c\ and exponent are com- 
putable functions of the ultraviolet scaling dimension of the fermion field. The AdS 
theory only considers a circular Fermi surface, and for this geometry (after appro- 
priate rescaling) 

v F (\k\-k F ) = v F (\q + k \-k F )^ q x + #, (65) 

now Eq. ( |64l is seen to be strikingly similar to Eq. ( |63l . Ref. l70l also argued that 
Eq. d64l i was a generic property of the near horizon geometry of the Reissner- 
Nordstrom black hole: the geometry changes from AdS4 near the boundary to 
AdS? x R 2 near the black hole horizon. 

It is interesting to compare the structure of the critical theory in the AdS/CFT 
framework to that found in the subsections above for the Ising-nematic transition in 
a metal. The latter was described by an infinite set of 2+1 dimensional field theories 
labeled by pairs of momenta on a one-dimensional Fermi surface i.e. a 5 • j2q set 
of 2+1 dimensional field theories. In the low-energy limit, the AdS/CFT approach 
yields [70] a AdS2 x R 2 geometry: this can be interpreted as an infinite set of chiral 
1+1 dimensional theories labeled by a R 2 set of two-dimensional momenta k. It is 
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notable, and perhaps significant, that both approaches have an emergent dimension 
not found in the underlying degrees of freedom. The Ising nematic theory began 
with a 2+1 dimensional Hamiltonian 5? c + + S? c § m Eqs. (14 1147148 1 ), and ended 
up with a S 1 /Zj set of 2+1 dimensional field theories. In AdS/CFT, there is the emer- 
gent radial direction representing energy scale. These emergent dimensions imply 
redundant descriptions, and require associated consistency conditions: we explored 
such consistency conditions in Section [{01 while in AdS/CFT the consistency con- 
ditions are Einstein's equations representing the renormalization group flow under 
changes of energy scale. It would be interesting to see if fluctuations about the clas- 
sical gravity theory yield corrections to the AdS2 x R 2 geometry which clarify the 
connection to our Ising-nematic theory. 
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